Angular momentum conservation for dynamical black holes 
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Angular momentum can be defined by rearranging the Komar surface integral in terms of a twist 
form, encoding the twisting around of space-time due to a rotating mass, and an axial vector. If 
the axial vector is a coordinate vector and has vanishing transverse divergence, it can be uniquely 
specified under certain generic conditions. Along a trapping horizon, a conservation law expresses 
the rate of change of angular momentum of a general black hole in terms of angular momentum 
densities of matter and gravitational radiation. This identifies the transverse-normal block of an 
effective gravitational-radiation energy tensor, whose normal-normal block was recently identified 
in a corresponding energy conservation law. Angular momentum and energy are dual respectively 
to the axial vector and a previously identified vector, the conservation equations taking the same 
form. Including charge conservation, the three conserved quantities yield definitions of an effective 
energy, electric potential, angular velocity and surface gravity, satisfying a dynamical version of 
the so-called first law of black-hole mechanics. A corresponding zeroth law holds for null trapping 
horizons, resolving an ambiguity in taking the null limit. 

PACS numbers; 04.70.Bw, 04.30.Nk 



I. INTRODUCTION 

The theory of black holes appears finally to be reach- 
ing a stage of maturity in which it can be applied in the 
most interesting, distorted, dynamic situations, with ap- 
propriate definitions of relevant physical quantities. This 
article mainly concerns angular momentum and its con- 
servation, which is the last major piece of what seems to 
be an essentially complete new paradigm for black holes. 
It therefore seems timely to review below the key ideas 
and results of what might be called the heroic, classical 
and modern eras. 

The first solution which would nowadays be called a 
black hole was discovered by Schwarzschild [1] within a 
few weeks of the final formulation of the field equations of 
General Relativity by Einstein 2], as the external grav- 
itational field of a point with mass M. Charge Q was 
soon added by Reissner [2i] and Nordstrom [3], but even 
the Schwarzschild solution was not properly understood 
for decades. Schwarzschild described the mass point as 
located at what is now understood as the horizon, despite 
its non-zero area A. Einstein & Rosen Q realized that 
the spatial geometry had a wormhole structure, extend- 
ing through a minimal surface. Oppenheimer & Snyder 
constructed a model of stellar collapse in which the 
star collapses through the horizon. Finally Kruskal Jl, 
as reported in a paper actually written by Wheeler Q, 
described the entire space-time, whence it became clear 
that there was a trapped region inside the horizon, from 
which nothing could escape to the exterior. Angular mo- 
mentum J was added by Kerr and, including Q, by 
Newman et al. [l^l- Uniqueness theorems identify these 
as the only black holes which are stationary, asymptoti- 
cally flat, electrovac solutions. 

Wheeler is credited with coining the term "black 
hole" and Penrose [l^l with defining event horizons, 
which became the accepted definition of black holes. 



Hawking [T^ showed that the area of an event horizon 
was non-decreasing, A' > 0. The result became known 
as a "second law" , due to inaccurate analogies with the 
laws of thermodynamics and the results summarized by 
Bardeen et al. [ij]: a "zeroth law" that surface gravity 
K is constant on a stationary black hole; a "first law" 
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for perturbations of stationary black holes, where Q is the 
angular speed and $ the electric potential of the horizon, 
and the ADM energy E measures the total mass of the 
space-time; and a "third law" that k 7^ by positive- 
energy perturbations of stationary black holes. This sum- 
marizes the classical theory of black holes as described in 
textbooks [ll,[li,[l3. 

The last results above are perhaps best described as 
black-hole statics, being properties of stationary black 
holes, specifically of Killing horizons, rather than of gen- 
eral event horizons. While adequate in some astrophysi- 
cal situations, this classical theory is inapplicable to gen- 
eral dynamical processes, for instance black-hole forma- 
tion, rapid evolution and binary mergers. A theory of 
black-hole dynamics is needed, with generalizations of all 
the above-mentioned quantities. Event horizons are not 
an appropriate platform, since they cannot be located by 
mortals, let alone admit physical measurements. A more 
practical way to locate a black hole is by a marginal sur- 
face, an extremal surface of a null hypersurface, where 
light rays are momentarily caught by the gravitational 
field. Marginal surfaces are used extensively in numeri- 
cal simulations, where they have historically been called 
apparent horizons^though the textbook definition of the 
latter is different Here a hypersurface foliated 

by marginal surfaces will be called a trapping horizon. 

A systematic treatment of trapping horizons [l^ dis- 
tinguished four sub-classes, called future or past, outer 
or inner trapping horizons, with the future outer type 
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proposed as the practical location of a black hole. Such 
a horizon was shown to have several expected proper- 
ties of a black hole, assuming the Einstein equation and 
positive-energy conditions: the horizon is achronal, being 
null in a special case of quasi-stationarity, but otherwise 
being spatial; the area form *1 of the marginal surfaces 
is constant in the null case and increasing in the spatial 
case; and the marginal surfaces have spherical topology, 
if compact. The area law implies that the area A — 
of the marginal surfaces S is non-decreasing. 



L^A > 
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where L denotes the Lie derivative and ^ the generating 
vector of the marginal surfaces. So a black hole grows if 
something falls into it, otherwise staying the same size. 

Comprehensive treatments were subsequently given for 
spherical symmetry [l^, |2l|, cylindrical s y mm etry 
[2^ and a quasi-spherical approximation (23l.[24[25l[2g|. 
In each case, definitions were found for all the non-zero 
physical quantities mentioned above, providing proto- 
types of all except J and O. In addition, an effective 
energy tensor Q for gravitational radiation was found, en- 
tering equations additively to the matter energy tensor T. 
The Einstein equations were decomposed into forms with 
manifest physical meaning, such as a quasi-Newtonian 
gravitational law, a wave equation for the gravitational 
radiation, and an energy conservation law which can be 
written in the form 



(3) 



where r in the normal dual of ^ and fc is a certain vec- 
tor, playing the role of a Killing vector, which is null on 
the horizon. Such equations actually hold not just on a 
trapping horizon, but anywhere in the space-time, energy 
conservation reducing at null infinity to the Bondi energy 
equation. 

Contemporaneousl y, A shtekar et al. and others [27l[28l. 
[2l,[33,[3ll,l3l[3l[3lll developed a theory of null trap- 
ping horizons with various additional conditions, under 
the names non-rotating isolated horizons, non-expanding 
horizons, weakly isolated horizons, rigidly rotating hori- 
zons and (strongly) isolated horizons. Each is intended 
to capture the idea that the black hole is quasi-stationary 
in some sense. They gave definitions of all the relevant 
physical quantities and derived a generalized version of 
the so-called first law. 

Subsequently, Ashtekar & Krishnan [3^,[33,|3^ studied 
future spatial trapping horizons under the name dynam- 
ical horizons, giving classes of definitions of energy and 
angular momentum, deriving corresponding flux equa- 
tions and obtaining a version of the so-called first law 
for (5 = 0. However, the "3-1-1" formalism used to de- 
scribe spatial trapping horizons breaks down when the 
horizon becomes null, so that the isolated-horizon and 
dynamical-horizon frameworks were essentially distinct. 
Some connections were drawn, particularly for slowly 



evolving horizons by Booth & Fairhurst [39|, |40|, |4l| . Re- 
cently, Andersson et al. [i^l showed that a stable trap- 
ping horizon is, on any one marginal surface, either spa- 
tial or null everywhere, so that transitions between the 
two types happen simultaneously on a marginal surface. 
They and Ashtekar & Galloway [43| also obtained some 
existence and uniqueness results for trapping horizons. 

A unified framework for any trapping horizon is pro- 
vided by a dual-null formalism 4Jj_ 45[, which was used 
throughout the earlier studies [iiniliO, HH, [H, [H, [13, 
[25I , [26j . The energy flux equation was then cast in a 
surface-integral form where the null limit could be taken 
(46l 1131 . Moreover, it was cast in the form of a conser- 
vation law ([3]), by identifying an effective energy tensor 
Q. The mass M, which might take any value on a given 
S by choice of scaling of fc, was chosen to be the irre- 
ducible mass or Hawking mass [i^ for consistency with 
the earlier studies. This corresponds to the simplest gen- 
eral definition of fc, such that it becomes a unit vector 
for round spheres near infinity. 

The main task here is to make similar refinements for 
angular momentum, as briefly described in shorter arti- 
cles 2^ [s^]. In particular, one desires not just a flux 
equation but a conservation law of the form 



LfJ 
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Here ijj should be an axial vector in some sense, playing 
the role of an axial Killing vector, with J being the angu- 
lar momentum about that axis. It turns out that natural 
restrictions on ip allow it to be uniquely specifled un- 
der certain generic conditions. The angular momentum, 
initially a functional J[V'], is obtained directly from the 
Komar integral [5l| in terms of a 1-form iv known as the 
twist [4^, which reduces to the 1-form used for dynami- 
cal horizons [H, [13, [ll] . It encodes the rotational frame- 
dragging predicted in the Lense-Thirring effect, thereby 
giving a precise meaning to the twisting around of space- 
time due to a rotating mass. 

The null limit is more subtle for angular momentum 
than for energy, where the irreducible mass is uniquely 
deflned for a null trapping horizon, energy conservation 
([3]) reducing correctly to L^M — 0. The dual- null folia- 
tion becomes non-unique for a null hypersurface, with 
uj becoming non-unique. This is reflected in the fact 
that different 1-forms were used for isolated horizons 
H [H [U, [13, [H, [13, [H and dynamical horizons 
. Neither 1-form is necessarily preserved along 
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a null trapping horizon, but a certain linear combination 
is so preserved. However, all three 1-forms coincide if the 
non-uniqueness of the dual-null foliation for a null hy- 
persurface is fixed in a certain way. Then J[ip] becomes 
unique for a null trapping horizon and conservation Q 
reduces as desired to L^J = 0. Thus a consistent treat- 
ment of angular momentum naturally resolves the issue 
of the degeneracy of the null limit. It follows that a 
black hole cannot change its angular momentum without 
increasing its area. 
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FIG. 1: A non-null hypersurface H foliated by spatial surfaces 
S, with generating vector ^ and its normal dual t = . \i H 
becomes null, ^ and r coincide. 

The article is organized as follows. Section II summa- 
rizes the underlying geometry. Section III reviews trap- 
ping horizons and conservation of energy. Section IV 
derives angular momentum from the Komar integral and 
shows how restrictions on the axial vector can be used to 
construct a unique definition. Section V derives and dis- 
cusses the conservation law. Section VI includes charge 
conservation and discusses local versus quasi-local con- 
servation. Section VII describes the state space, defining 
the remaining physical quantities and deriving a dynami- 
cal version of the so-called first law. Section VIII consid- 
ers null trapping horizons, deriving a zeroth law. Appen- 
dices concern (A) weak fields, (B) normal fundamental 
forms and (C) a Kerr example. The above discussion 
serves as a summary. 



II. GEOMETRY 

General Relativity will be assumed, with space-time 
metric g. The geometrical object of interest is a one- 
parameter family {S} of spatial surfaces S, locally gener- 
ating a foliated hypersurface H. Labelling the surfaces by 
a coordinate x, they are generated by a vector ^ = d/dx, 
which can be taken to be normal to the surfaces, = 0, 
where _L denotes projection onto S. A Hodge duality 
operation on normal vectors r/, 1.7] ~ 0, yields a dual 
normal vector rj* satisfying 

-L??* = 0, g{Ti*,-n) = 0, g{v*,V*) = ^giv, v)- (5) 
In particular, 

r = r (6) 

is normal to H, with the same scaling (Fig[T]). The co- 
ordinate freedom here is just x i— > x(x) and choice of 
transverse coordinates on 5', under which all the key for- 
mulas will be invariant. The generating vector ^ may 
have any causal character, at each point. For instance, 
a future outer trapping horizon is spatial while growing, 
becomes null when quasi-stationary, and would become 
temporal if shrinking during evaporation p^ . |46| . |47| . 

A dual-null formalism [44, 45] describes two families of 
null hypersurfaces E±, intersecting in a two-parameter 
family of spatial surfaces, including the desired one- 
parameter family. Some merits of the dual-null approach, 
apart from comparative ease of calculation, are that it is 
adapted both to marginal surfaces, defined as extremal 



surfaces of null hypersurfaces, and to radiation propaga- 
tion, which makes it easier to identify terms arising due 
to gravitational radiation [i^, [13, Ho] ■ Relevant aspects 
of the formalism are summarized as follows. 

Labelling by coordinates x^ which increase to 
the future, one may take transverse coordinates x"^ on 
5, which for a sphere would normally be angular co- 
ordinates x"^ = ("djip). Writing space-time coordinates 
x" = {x'^,x~,x°') indicates how one may use Greek let- 
ters (q;,/3, . . .) for space-time indices and corresponding 
Latin letters (a, b, . . .) for transverse indices. The coordi- 
nate basis vectors are da — d/dx"" and the dual 1-forms 
are dx", satisfying df3{dx°') = J^g. Coordinate vectors 
commute, [da, dp] ~ 0, where the brackets denote the 
Lie bracket or commutator. Two coordinate vectors have 
a special role, the evolution vectors d± — d/dx^ which 
generate the dynamics, spanning an integrable evolution 
space. The corresponding normal 1-forms dx^ are null 
by assumption: 

g-\dx^,dx^) ^0. (7) 

The relative normalization of the null normals may be 
encoded in a function / defined by 

^ -g-\dx+,dx-) (8) 

where the metric sign convention is that spatial metrics 
are positive definite. The transverse metric, or the in- 
duced metric on S, is found to be 

h = g + 2e^f dx+ (g) dx~ (9) 

where (gi denotes the symmetric tensor product. There 
are two shift vectors 

s± = ±5± (10) 

where _L is generalized to indicate projection by h. The 
null normal vectors 

l±^d±-s± = -e-^'g''\dx^) (11) 

are future-null and satisfy 

g{l±J±) = 0, g{l+,l-)^-e-f, 

l±{dx^) = 1, /±(da;=F) = 0, ±l± = 0. (12) 

The metric takes the form 

g = habidx'^ + s'^dx'^ + s°Ldx~) <^ 

{dx'' + sXdx+ + stdx-) - 2e-fdx+ ® da;~.(13) 

Then (/i, /, s±) are configuration fields and the indepen- 
dent momentum fields are found to be linear combina- 
tions of the following transverse tensors: 



6'± = *L±*1 (14) 

cr± = ±L±h-e±h (15) 

iy± = L±f (16) 

u = ie^M[L,M) (17) 
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where * is the Hodge operator of h and L± is shorthand 
for the Lie derivative along /-|-. Then the functions 9± 
are the null expansions, the traceless bilinear forms <j± 
are the null shears, the 1-form uj is the twist, measuring 
the lack of integrability of the normal space, and the 
functions iy± are the inaffinities, measuring the failure of 
the null normals to be affine. The fields {9±,(T±,h'±,uj) 
encode the extrinsic curvature of the dual-null foliation. 
These extrinsic fields are unique up to interchange ± i— > 
=F and diffeomorphisms i-^ x^{x^) which relabel the 
null hypcrsurfaces. Further description of the geometry 
was given recently [47|. 

As described, the dual-null formalism is manifestly co- 
variant on S, with transverse indices not explicitly de- 
noted, while ± indices indicate the chosen normal ba- 
sis (45l . . Conversely, one can use a formalism which 
is manifestly covariant on the normal space, with trans- 
verse but not normal indices explicitly denoted [sS] . Both 
types of formalism can seem obscure to the uninitiated, 
so indices will be explicitly denoted in longer formulas in 
this article, nevertheless being omitted where the mean- 
ing should be clear. Capital Latin letters (A, B, . . .) will 
be used for normal indices, when not denoted by ± in 
the dual-null basis. Then the configuration fields are 
{hab,f,SA''), the momentum fields are {9A,crAbc,i^A,i^a) 
and the derivative operators are (J-La, Da), where D is 
the covariant derivative operator of h. 

Since the normal space is not integrable unless uj = 0, 
it generally does not admit a coordinate basis. However, 
one may still take dx^ as basis 1-forms, in which case the 
dual basis vectors are l±, as follows from ([T^ . implying 
lA{dx^) = 5^. In this basis, the normal metric 

l = 9-h (18) 

has components which follow from ([9]) as 

Jab = —e^^ {dx^dxg + dx^dx^) (19) 

and its inverse has components 

= -ef (l^l^^ + l^lf) (20) 

which can be used to lower and raise normal indices. Also 
useful is the binormal 

£AB = e^-^{dx\dxg — dx^dx^) (21) 

or its inverse 

^^f{lAiB_lAiBy ^22) 

The mixed form 

e^B = lidx% - l^dx-g (23) 

has components = ±1, e^zp = 0, so can be used 
to express the duality operation ([5|) on normal vectors, 
extended to the dual-null foliation, as 



The dual-null Hamilton equations and integrability 
conditions for vacuum Einstein gravity were derived pre- 
viously [i^, with matter terms added subsequently [25j . 
The components of the field equations which are rele- 
vant to angular momentum turn out to be the twisting 
equations 

± \h^''Dda±ac T 87rT,± (25) 

where Ta± = hlT^pl^ is the transverse-normal projec- 
tion of the energy tensor T, and units are such that New- 
ton's gravitational constant is unity. The corresponding 
all-index version can be written using the binormal as 

LLBUJa = -9BtOa + \e''B{DaVE- DjE~9EDaf 

+ h'^^DdOEac - IGl^TaE). (26) 

III. TRAPPING HORIZONS AND 
CONSERVATION OF ENERGY 

Returning to a general foliated hypersurface iJ, a nor- 
mal vector r/ has components ry^ along Z±, so that rj = 
ri^l+ + rj~l-, and its normal dual is rj* = rj~^l+ — rj^l-. 
In particular, the generating vector is 

(27) 

and its dual is 

T^ei+-ri- (28) 

Since the horizon is given parametrically by functions 
x'^{x), the components = dx'^/dx are independent of 
transverse coordinates: 

= 0. (29) 
It is also useful to introduce the expansion 

9, = L.logi^l) = 9av'^ (30) 

along a normal vector 77, particularly the expansion 9^ 
along the generating vector. 

A trapping horizon [H, IH, \^ is a hypersurface H 
foliated by marginal surfaces, where S is marginal if one 
of the null expansions, 9^ or vanishes everywhere on 
S. Then S is an extremal surface of the null hypersurface 
S+ or 

The recently derived energy conservation law fi^, |43| 
will be stated here for later comparison, modifying some 
notation. Assuming compact S henceforth, the trans- 
verse surfaces have area 

j> *l (31) 

and the area radius 



(24) 



R = V^/47r 



(32) 



5 



is often more convenient. The Hawking mass [4a | 

can be used as a measure of the active gravitational mass 
on a transverse surface. Assuming the null energy con- 
dition, M is the irreducible mass R/2 of a future outer 
trapping horizon, L^AI > 0, by the area law On 
a stationary black-hole horizon, M reduces to the usual 
definition of irreducible mass for a Kerr-Newman black 
hole, namely the mass which must remain even if rota- 
tional or electrical energy is extracted. It is generally not 
the ADM energy, but an effective energy E is defined in 
§VII which does recover the ADM energy in this case. 
Equality on a trapping horizon will be denoted by 
e.g. R^2M. 

Mass or energy has a certain duality with time, e.g. 
there is a standard formula for energy if a stationary 
Killing vector exists. For a general compact surface, the 
simplest definition of such a vector which applies cor- 
rectly for a Schwarzschild black hole is [i^, l47j 

k^{g-\dR))* (34) 

or — e^^LsR- This vector actually was found to 
be the appropriate dual of M, in the sense of conserva- 
tion laws for trapping horizons [ 46l , [47| and for uniformly 
expanding flows [52, 53]. In cither case, the energy con- 
servation law can be written as 

L^M = j^*{TAB + QAB)k^T^ (35) 

where O is an effective energy tensor for gravitational 
radiation. This determines only the normal-normal com- 
ponents of O, as 

e±± = lk±||V32^ (36) 
e±^ = e-^\uoT\Df\^l&^ (37) 

where |CP = h'^^CaCb andilcrlp = h'"'h'"^aab<Tcd. Further 
discussion is referred to 0, 153, [H, . 

IV. ANGULAR MOMENTUM 

The standard definition of angular momentum for an 
axial Killing vector ip and at spatial infinity is the Komar 
integral [5l| 

=-Y^£«a/3V"7A''. (38) 

Now consider ip tohe a, general transverse vector, J-ip = ip 
(FiglJ). Since taptp^ = 0, the Komar integral can be 
rewritten via as 

JVl^] (39) 
f^TT JS 




FIG. 2; A transverse vector ip. 

where lu is the twist (flT)) . Since the twist encodes the 
non-integrability of the normal space, it provides a geo- 
metrical measure of rotational frame-dragging. It is an 
invariant of a dual-null foliation and therefore of a non- 
null foliated hypersurface H, so the twist expression for 
J[ip] is also an invariant. Appendix A shows that J[^] re- 
covers the standard definition of angular momentum for 
a weak- field metric (IQi] , with the twist being directly re- 
lated to the precessional angular velocity of a gyroscope 
due to the Lense-Thirring effect. Thus the twist does in- 
deed encode the twisting around of space-time caused by 
a rotating mass. 

There are several definitions of angular momentum 
which are similar surface integrals of an axial vector 
contracted with a 1-form [3^, l3a. [s^ l , the situation be- 
ing clarified by Gourgoulhon [55| and in Appendix B. 
Ashtekar & Krishnan [3y] gave a definition for a dynam- 
ical horizon which involves a 1-form coinciding with uj. 
Brown & York [5^] gave a definition which was stated 
only for an axial Killing vector -0, involving a 1-form 
which is generally inequivalent to uj, but can be made to 
coincide if adapted to a trapping horizon. Ashtekar et 
al. [sO, HH gave a definition for a type II (rigidly rotat- 
ing) isolated horizon, using a 1-form which is generally 
inequivalent to uj. However, it can be made to coincide 
with UJ if the dual-null foliation is fixed in a natural way, 
as described in the penultimate section. 

The above properties suggest ((39)) as a general quasi- 
local definition of angular momentum. However, if the 
transverse vector tp does not have properties expected 
of an axial vector, the physical interpretation as angu- 
lar momentum is questionable. For instance, it would be 
natural to expect an axial vector to have integral curves 
which form a smooth foliation of circles, apart from two 
poles, assuming spherical topology for S. In the follow- 
ing, two conditions on ip with various motivations will be 
considered, which, taken together, determine %p uniquely 
in a certain generic situation. These conditions then yield 
a conservation law with the desired form Q , as described 
in the next section. 

Ashtekar & Krishnan 37] proposed that ip has vanish- 
ing transverse divergence: 

2?aV'° = 0. (40) 

This condition holds if ip is an axial Killing vector, and 
can be understood as a weaker condition, equivalent to 
Ip generating a symmetry of the area form rather than 
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FIG. 3: Curves 7 of constant expansion d^. 

of the whole metric, since L^{*1) ~ ^Daip"". Aherna- 
tively, assuming that the integral curves of ip are closed, 
it can always be satisfied by choice of scaling of V^, as 
discussed by Booth & Fairhurst [13] ■ The original moti- 
vation was that the different 1-forms used for dynamical 
and isolated horizons, denoted here by lu and lu + \Df, 
will then give the same angular momentum J [?/'], by the 
Gauss divergence theorem. 

Spherical topology will be assumed henceforth for S, 
which follows from the topology law [l^ for outer trap- 
ping horizons, assuming the dominant energy condition. 
If there exist angular coordinates (1?, Lp) on S with ■0 = 
d/dip, completing coordinates {x,'d,ip) on H, then since 
coordinate vectors commute, 

L^V = 0. (41) 

This condition was proposed as a natural way to prop- 
agate -0 along H by Gourgoulhon [s^. Now there is a 
commutator identity (45| 

L^{Dar ) ~ DaiL^i'T = ^DaO^. (42) 

Therefore assuming both conditions (|l(Il) - (|¥T|l forces 

^''Da0^ = 0. (43) 

This is automatically satisfied if D9^ = 0, as in spheri- 
cal symmetry or along a null trapping horizon. However, 
generically one expects D9^ ^ almost everywhere. It 
must vanish somewhere on a sphere, by the hairy ball 
theorem, but the simplest generic situation is that there 
are curves 7 of constant 9^ which form a smooth folia- 
tion of circles, covering the surface except for two poles 
(FiglSI). Assuming so, since V is tangent to 7, one can 
find a unique 7/;, up to sign, in terms of the unit tangent 
vector ^ and arc length ds along 7: 

■ip^^p (b ds/2TT (44) 

where the scaling ensures that the axial coordinate is 
identified at and 2tt. Then the angular momentum 
becomes unique up to sign, J[ip] — J, the sign being nat- 
urally fixed by J > and continuity of tp, corresponding 
to a choice of orientation. 

For an axisymmetric space-time with axial Killing vec- 
tor ijj, (|40p holds, while (|4ip holds if ^ respects the sym- 
metry, L^^ = 0, so the above construction, if unique as 
assumed, yields the correct axial vector. In particular. 



the construction does work for a Kerr space-time, as de- 
scribed in Appendix C. 

To summarize this section, the definition ([55)1 of an- 
gular momentum can be made generically unique if the 
axial vector is a coordinate vector, (|4T|) . and generates a 
symmetry of the area form, (|40p . The construction can 
be applied in any situation where D9(^ ^ almost ev- 
erywhere, though the physical interpretation as angular 
momentum seems to be safest in the case of two poles, 
which locate the axis of rotation. Then J is proposed to 
measure the angular momentum about that axis. 

V. CONSERVATION OF ANGULAR 
MOMENTUM 

The main result of this paper is that 

L^J^- j^^iTaB + QaBWr'' (45) 

holds along a trapping horizon under the conditions (HDl) - 
(|4T|l . where 

QaB - -^h^'-Ddasac (46) 

is thereby determined to be the transverse-normal block 
of the effective energy tensor for gravitational radiation. 
It can be shown by differentiating the angular momentum 
([39)1 using gH) to give 

L^J^^j^*{9^^''Ua+rLiU0a) (47) 

then expanding ^ by (P?]) and using the twisting equa- 
tions ([25]) to express L^uj. The term in 9^ = £^'^9++^^9^ 
cancels with the first term from (pS)) . while the D gradi- 
ents may all be removed as total divergences due to (^5]) , 
(gni) and the fact that (gS]) reduces to tp''Da9^ = on a 
trapping horizon with 9± ^ 0. This leaves just 

-n^h^'^DdCT.ac - 87rT,_))(48) 

which is an expanded form of noting and 

thereby identifying (|46|) . 

Apart from the inclusion of O, the conservation law 
(j45p is the standard surface-integral form of conservation 
of angular momentum, were ip an axial Killing vector. 
It thereby describes the increase or decrease of angular 
momentum of a black hole due to infall of co-rotating or 
counter-rotating matter respectively. The corresponding 
volume-integral form for a spatial horizon H, expressing 
the change [J]dH in J between two marginal surfaces, 
follows as 

[J]0H = - / HTaB + QaBWr'' (49) 
JH 
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where f = t/ yjg{^, is the unit normal vector and *1 = 
g{i, A dec is the proper volume element. Although 
more familiar, as for the energy conservation law (46l.[?n|. 
this form becomes degenerate as the horizon becomes 
null, since *1 — > while f ceases to exist. Since this is 
a physically important limit, where a black hole is not 
growing, the surface-integral form ()45|) is preferred. 

The null shears a±hc have previously been identified in 
the energy conservation law (|35p as encoding the ingo- 
ing and outgoing transverse gravitational radiation, via 
the energy densities p6p . which agree with expressions in 
other limits, such as the Bondi energy density at null in- 
finity and the Isaacson energy density of high-frequency 
linearized gravitational waves [i^, [13] . So the expression 
(I46p implies that gravitational radiation with a trans- 
versely differential waveform will generally possess angu- 
lar momentum density. One can see corresponding terms 
in the linearized approximation fl6|, but they are set to 
zero by the transverse-traceless gauge conditions, which 
are "transverse" in a different sense to that used here. In 
any case, the conservation law shows that a black hole 
can spin up or spin down even in vacuum, at a rate re- 
lated to ingoing and outgoing gravitational radiation. 

The identification of the transverse-normal block ((46l) 
of appears to be new. Previous versions of an- 
gular momentum flux laws for dynamical black holes 
[Hi [13, m, m, SO, SH, ill contain different terms, 
which are not in energy-tensor form, i.e. some ten- 
sor contracted with ■0 and r. They can be recovered 
by removing a transverse divergence from OqsV'"''"^, 
yielding afDd^pa/l&T^ = <'^i^/iad/327r, where af = 
T^h°''^h"^aBce encodes the shear along r. Such terms 
have been described by analogy with viscosity (55 . 56] . 

The conservation laws ([55]) and (|45p take a similar 
form, expressing rate of change of mass M and angular 
momentum J as surface integrals of densities of energy 
and angular momentum, with respect to preferred vec- 
tors k and ip which play the role of stationary and axial 
Killing vectors, even if there are no symmetries. Of the 
ten conservation laws in flat-space physics, they are the 
two independent laws expected for an astrophysical black 
hole, which defines its own spin axis and centre-of-mass 
frame, in which its momentum vanishes. 



VI. QUASI-LOCAL CONSERVATION LAWS 

For an electromagnetic field with charge-current den- 
sity vector j, the total electric charge Q in a region H of 
a spatial hypersurface is defined as 



[Q]aH = - / *g{i,T) 

JH 



(50) 



The surface-integral form of conservation of charge fol- 
lows by the same arguments relating (|45p and 



As before, this is more general, since H may have any 
signature. The conservation laws for energy (|35p and 
angular momentum ()45p evidently take the same form 



L^M = -6 *g(J, r), L^J = - f *g{j, t) (52) 
Js Js 

by identifying current vectors 

f = -kA{T^^ + j^ = Va(T'^'' + e"^). (53) 

The local differential form of charge conservation. 



(54) 



where V is the covariant derivative of g, notably does 
not hold for J or j in general. A weaker property holds, 
obtained as follows. First note that in any of the three 
conservation laws ([^5)1 and ([5T|) . ^ and r may be in- 
terchanged. Thus there are really two independent laws 
in each case. This can be understood from special rel- 
ativity: if ^ were causal, one would interpret them as 
expressing rate of change of energy, angular momentum 
or charge as, respectively, power, torque or current; while 
if ^ were spatial, one would normally convert to volume- 
integral form and regard them as defining the energy, 
angular momentum or charge in a region. One can make 
either interpretation for a black hole, since H would be 
generically spatial, but the marginal surfaces S locate the 
black hole in a family of time slices. 

Given two independent equations in the normal space, 
it follows that 



e-^siAM = - * *1 A Jb, e'^BLAJ = -f*l/\jB- 
Js Js 

(55) 

Expressed in terms of the curl and divergence of the nor- 
mal space, 



curl M • 
whereas 



*1AJ, curlJ = - a *1 A J (56) 
s Js 



V^f - *div(*l A j) (57) 
for any normal vector ]. Then divcurl — yields 

(58) 



(51) 



*V„r = f *^af = 0. 



This subtly confirms the view that energy and angular 
momentum in General Relativity cannot be localized 1™ , 
but might be quasi- localized, as surface integrals 57|. 
The corresponding conservation laws have indeed been 
obtained in surface-integral but not local form. 



VII. STATE SPACE 

There are now three conserved quantities {M,J,Q), 
forming a state space for dynamical black holes. Fol- 
lowing various authors [ll|, [H, HO, 113, [H, HI], related 
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quantities may then be defined by formulas satisfied by 
Kerr-Newman black holes, specifically those for the ADM 
energy 



E 



V((2A/)2+Q2)2 + (2J)2 



the surface gravity 



(2M)4 - (2J)2 - Q4 



2(2M)3v/((2Af)2 + Q2)2 + (2j)2 



the angular speed 



n = 



J 



M^((2Af )2 + Q2)2 + (2 j)2 



and the electric potential 



((2M)2 + Q2)Q 



2Af^((2M)2 + g2)2 + (2J)2 



(59) 



(60) 



(61) 



(62) 



It would be preferable to have independently motivated 
definitions of these quantities, but so far this has been 
done only in spherical symmetry [l^, [l^, [lH , where there 
are natural definitions of E, k and ^ = Q/R which can 
be applied anywhere in the space-time, coinciding with 
the above expressions on the outer horizons of a Reissner- 
Nordstrom black hole. 

In the dynamical context, E > M is generally not the 
ADM energy, since there may be matter or gravitational 
radiation outside the black hole. Rather, it can be in- 
terpreted as the effective energy of the black hole, as 
follows. Defining the moment of inertia / by the usual 
formula J = lil yields 



/ ^ M V((2M)2 -hQ2)2 + (2j)2 5. . (63) 

Expanding for J < and Q < M yields 

£; w M + + \Q^/R. (64) 

The second and third terms are standard expressions for 
rotational kinetic energy and electrostatic energy. Thus 
the irreducible mass M plays the role of a rest mass, with 
E including contributions from rotational and electrical 
energy. 

Returning to the general case, the above definitions 
satisfy the state-space formulas 

_dE_ I dE ^^dE ^^dE 

^ = ^^dA = ^idA-r ^=97' *=aQ- 

There follows a dynamic version of the so-called first law 
of black- hole mechanics Il4| : 



L(^E 



-L^A + QL^J + <i>L^Q. 



(66) 



As desired, the state-space perturbations in the classical 
law for Killing horizons or the version for isolated 



horizons [30|, |3l|, , have been replaced by the deriva- 
tives along the trapping horizon, thereby promoting it to 
a genuine dynamical law. 

The rate of change of effective energy can also be writ- 
ten in energy-tensor form, 

L^S = £ * {{To^p + ea/3)if"r'3 - ^j^t/^) (67) 



where 



K = AMnk - VliP 



(68) 



plays the role of the stationary Killing vector. Note that 
in the classical theory of stationary black holes, the state 
variables are usually taken to be (i?, J, Q), with the ir- 
reducible mass M defined as a dependent variable. The 
dynamical theory reveals that (Af, J, Q) are the more ba- 
sic variables, since they each satisfy a simple conservation 
law. Then the effective energy E is defined as a depen- 
dent variable and therefore satisfies the above conserva- 
tion law. This refiects a shift in emphasis from the clas- 
sical to the dynamical theory: the so-called first law is a 
dependent result, obtained from more fundamental con- 
servation laws for energy, angular momentum and charge. 



VIII. NULL TRAPPING HORIZONS AND 
ZEROTH LAW 

A zeroth law for trapping horizons follows from the 
above, if one defines local equilibrium by the absence of 
relevant fluxes: 



(69) 



Then (Af , J, Q) are constant on the horizon and so is 
K. In fact, these conditions do hold on a null trapping 
horizon under the dominant energy condition, as shown 
below. This treatment also turns out to be compatible 
with the definition of weakly isolated horizon [2^ [3l|, [3^ . 

Consider a null trapping horizon, assumed henceforth 
in this section to be given by 0+ =0. The null focussing 
equation yields T++ -|- | |(t_)_ | p/327r = 0, so the null energy 
condition, which implies r++ > 0, yields [iB] 



0, (J, 



0. 



(70) 



Thus the degenerate metric of the horizon is preserved 
along the generating vector. The dominant energy con- 
dition, which implies that the energy-momentum Pa — 



Tapl^ is causal, further yields [2£ 



+a 







(71) 



since P = —T^+dx^—T^ dx^ —T+adx'^ would otherwise 

be spatial. 

On a null hypersurface H , one can take the null coor- 
dinate to be the generating coordinate, x^ = x, meaning 
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that the shift vector vanishes, s+ = 0, so that ^ = t = l^. 
Since k = —e^L^Rl^ (|34p . one finds 

5(J,t) = e/i_i?(T+++e++), g{j,T) = r{Ta++ea+). 

(72) 

These fluxes vanish by the above results ([70 |) -([7T |l and 
the expressions (|36|) . (|46l) for components of 9. The other 
flux in vanishes due to the Maxwell equations f29j . 

For a null trapping horizon, the dual-null foliation is 
not unique, so the question arises whether there is a nat- 
ural way to fix it. An affirmative answer is given by 
noting that the above results also imply, via the twisting 
equations (p5)). 

L+(u;-\Df)^G. (73) 

This restriction on the dual-null geometry is suggestive 
of a proto-conservation law for angular momentum. Now 
the only normal fundamental form intrinsic to a null hy- 

persurface is Q ^ = ^Df + uj oi A ppe ndix B, which was 

therefore used by Ashtekar et al. [23| to define angular 
momentum for an isolated horizon. However, it is the 
other null normal fundamental form C_| — ^Df — ui, de- 
pending on the dual-null foliation, which is preserved as 
above. Given that the general definition (|39p of angular 
momentum involves u, it seems best to fix the unwanted 
freedom by 

Df - 0. (74) 

Recalling the definition ([8]) or (fT2)) of /, this condition 
fixes the normalization of the extrinsic null normal Z_ 
with respect to the intrinsic null normal Z+, which is al- 
ways possible on a null hypersurface H. In fact, it is 
common simply to fix / = 0. A similar normalization is 
also used in the context of null infinity. 

Then the definition p9p of angular momentum be- 
comes unambiguous on a null trapping horizon, coincides 
with the definition for isolated horizons [s^, [3l|, [s^ , and 
is preserved along the horizon, assuming only that ip is a 
coordinate vector field (|¥T|l : 

L^J^O. (75) 

Since the area law (l^ shows that A is increasing unless 
H is null everywhere on a given S, this answers, in the 
negative, a simply stated physical question: whether a 
black hole can change its angular momentum without 
increasing its area. 

The above reasoning has largely recovered the notion 
of a weakly isolated horizon introduced by Ashtekar et al. 
[2^ . except that the scaling freedom in Ll. has not been 
fixed. In more detail, Ashtekar et al. [29[ introduced a 
1-form which will here be denoted by w, defined by 

VJI = vjJI (76) 

where V is the covariant derivative operator of H . The 
transverse and normal components are found as 

1.VJ = —u — jDf, I'^Wa = (77) 



Since w is an invariant of H and Ashtekar et al. (29| 
demanded 

L+vj ^ (78) 

in order to define a weakly isolated horizon. The trans- 
verse part agrees with the above results, which also imply 
Diy^ ~ 0, while the normal part further fixes the inaffin- 
ity to be constant on H. This fixes the scaling of /+ 
up to a constant multiple. Ashtekar et al. [2^ defined 
the surface gravity to be 

k ^ -v+ (79) 

which recovers the standard surface grav ity of a Killing 
horizon if is the null Killing vector [17| . Then the con- 
stancy of k can also be interpreted as a zeroth law. This 
still leaves non-zero k ambiguous up to a constant mul- 
tiple, not necessarily agreeing with the definition (j60p . 
which therefore fixes that freedom. 

The above considerations appear to have a closed a 
gap in the general paradigm, concerning how a growing 
black hole ceases to grow. It seems that the generically 
spatial trapping horizon simply becomes null. It is dif- 
ficult to find a practical formalism describing all cases 
without some degeneracy arising in the null case, but the 
dual-null formalism appears to be adequate; one fixes 
the additional freedom in the null case by ([74]) . In par- 
ticular, no additional conditions need be imposed on the 
horizon itself, as compared with the variety of definitions 
of isolated horizons [23, [H, [s^, [Ml, [H IH, [H, |3|. 
Numerical evidence that such horizons exist in practice 
has been given by Dreyer et al. [s^l, who looked for and 
found approximately null trapping horizons, under the 
name non-expanding horizons. 

Research supported by the National Natural Sci- 
ence Foundation of China under grants 10375081 and 
10473007 and by Shanghai Normal University under 
grant PL609. Thanks to Abhay Ashtekar, Ivan Booth, 
Eric Gourgoulhon and Badri Krishnan for discussions. 

APPENDIX A: TWIST AND WEAK FIELDS 

The twist may be calculated by first finding the shift 
vectors s±, due to the form [45] 

u: = \efh ([a+, s_] - s+] + s+]) . (Al) 

If it is more convenient to use an orthonormal basis 
{Zo, ll} of the normal space, 

_LZo = -L^i = = .g(/o, h), g{h, k) = -1 = -gihji) 

(A2) 

then 

uj = ^h{[lo,li]) (A3) 

follows by linear combinations from (I39p , or directly from 
the Komar integral (l38|) . If the basis is adapted to a 
coordinate basis via 

Ia = dA - SA, SA = -LOa (A4) 
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then 



cj = i/i So] - [do, si] + [so, si]) 



(A5) 



In either case, the first s tep is to find the shift vectors. 

The weak-field metric [l6| , in standard spherical polar 
coordinates {t, r, i9, (p) adapted to the axis of rotation, is 



2M 

r 
2M 

r 



4J . 2 

— sin "ddtdip 
r 

2 



{dr^ + r^{di9^ + sin^ Mip^)){A6) 



where, in this appendix only, M and J denote the mass 
and angular momentum as defined in this approximation, 
obtained by linearizing the metric and neglecting higher 
powers of l/r. The inverse metric is 



1 



4J 



dtd^ 




Taking the transverse surfaces S as those of constant 
{t, r), one can read off the nonzero component of the shift 
1-forms SAb as 



9tv 



2J 



sin^ 



St' 



2J 



(A8) 



Then the non-zero component of the twist is given by 
3J 3J 



■ sin^ i). 



The area form is 

*1 ^ sin-dd-d A dip 

so that 

*ujip ~ SJsin'^ A dip. 
Standard integrals yield 

-uju> ^ Stt J. 



(A9) 



(AlO) 



(All) 



Since 0Ja'4'°' = uj^p if ip — d/d(p, this agrees with the gen- 
eral definition (|39|) of angular momentum. 

A directly measurable quantity due to rotational 
frame-dragging is the precessional angular velocity 



f^LT ~ ^ (3(i • r)f - z) 



(A13) 



of a gyroscope due to the Lense-Thirring effect, where 
f is a unit vector in the direction of the gyroscope and 
z is a unit vector along the axis of rotation. Results of 
measurements of the effect due to the Earth by Gravity 



Probe B are expected soon. If the twist uj is formally 
converted to an angular velocity D, by 



to ^ Q X f 



then 



UJ 



sint? 



3J 



(A14) 



(A15) 



does have the direction and relativistic dimensions of an- 
gular velocity. Then 



n 



LT 



(A16) 



Curiously, this is a linear transformation of 17, by the 
same traceless tensor used in defining quadrupole mo- 
ments !l6j. The puckish role of the factor of 3 in the 
above calculations is also noteworthy. In any case, it 
confirms that the twist docs indeed encode the twisting 
around of space-time due to a rotating mass, in a directly 
measurable way. 

For completeness, the agreement of the weak-field mass 
with the Hawking mass can also be checked, as follows. 
One needs to keep track of an extra power of 1 /r in the 
area form 



*1 



1 



2Af 



sin ddd A dip. 



(A17) 



Then the expansion 1-form 9a has non-zero component 



o , , , 2M 
9, log 1 + 



Then 



and 



^ 1 



(A18, 



4 1 



s 

IGlT ( 1 



-) " 



sin ddi) A d(p 

(A20) 



(A12) Since r, this agrees with the Hawking mass (|33|) . 



APPENDIX B: NORMAL FUNDAMENTAL 
FORMS 

Various definitions of angular momentum [s^, Is^L [s^ l 
are similar to (|39p . with lo replaced by a 1-form which 
is, implicitly or explicitly, a normal fundamental form. 
To clarify the situation, normal fundamental forms are 
reviewed below, referring to previous treatments [45l [55| . 

Writing the twist (fT7|) explicitly in components, 



2u;o 



'i^Kp{ll^-,ll-llV-,f_). 



(Bl) 



11 



If l± are adapted to a coordinate basis via pT|) . then 
commutativity Vj^V^jx^ = allows it to be written as 



(B2) 



This is the difference of two normal fundamental forms 
Ct± with components 



7- 



(B3) 



They are the independent normal fundamental forms, 
since the corresponding C,±± vanish. Their sum is Df, 
since the normalization in (I12|l yields 



Then 



Do.f ^e^hHrLVpl+^ + llVfil-^). (B4) 



(B5) 



and the normal fundamental forms are thereby encoded 
in to and Df. 

For an orthonormal basis (jA2p . there is just one in- 
dependent normal fundamental form Coi — — CiOi given 
by 

Uc.=llhi^ph^, Cioc.=r,hyplo^ (B6) 

with corresponding C,qq, (.h vanishing. Under a boost 
transformation 

^0 coshp+Zi sinh/9, i— )■ sinhp+Zi coshp (B7) 



which preserves the orthonormal conditions (jA2p . the 
normal fundamental form is generally not invariant: 



Coi Cio - Dp, Cio ^ Cio + Dp. 



(B8) 



However, if the orthonormal basis is adapted to the dual- 
null basis by, e.g. V^l± = Iq + h, then 



Then 



Coi — ^Cio — ({ 



00 



-Cii = hD.f. 



(B12) 



These 1-forms nevertheless become degenerate when ^ 
becomes null. 

As shown explicitly by Gourgoulhon ^55*1, the 1-form 
used by Brown & York [5J] to define angular momentum 
is, in this notation, —Cio- Therefore it coincides with uj 
if the orthonormal basis is adapted to a trapping hori- 
zon, but generally not if it is adapted to a foliation of 
spatial hypersurfaces intersecting the trapping horizon 
in the marginal surfaces. The 1-form used by Ashtekar 
& Krishnan [36] to define angular momentum for dynam- 
ical horizons is also —Cio, this time explicitly adapted to 
the horizon by (jBlOp , therefore coinciding with to in this 
case. The 1-form used by Ashtekar et al. [s^, [U, [s^l to 

define angular momentum for isolated horizons is C hi 

which generally does not coincide with co. However, it 
does coincide if the gauge freedom is fixed by ((74)) . 



APPENDIX C: KERR EXAMPLE 



Consider a Kerr space-time in Boyer-Lindquist coor- 
dinates (i, r, (/j), with S given by constant {t,r) and 
C = d/dr. Then 



1 = J (r2 + a2)2 - Aa2 sin^ sin Mi!) A d(p (CI) 



where A 



2mr + a^, 



Coi — —Cio — 



(B9) 



For spatial H, the same result is obtained by adapting 
the orthonormal basis by choosing 



_ 2r(r2 -t- a2) - (r - m)a2 sin^ i) 

5 ^ 'n~, 2^2 A 2 • 2 Q y^^) 

(r^ -I- a^)^ — Aa^ sm tf 



(BIO) and 



The missing information in Df can be recovered by in- 
stead defining 



COIQ 

CoOct 



<,10a - 



(Bll) 



2a2(r2 + a2)(r^ - 3mr2 -|- a^r + ma^) sin i? cost? 

Due — r; dij. 

5 ((r2 + a2)2 - Aa2 sin^ ,?)2 

(C3) 

If a 7^ 0, D6^ is non-zero except at the poles, equator and 
isolated values of r, so the construction yields a unique 
continuous ip, ip = d/dip. 
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